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1 Introduction 



In recent times there has been an enormous amount of work devoted to the 
understanding of random surfaces and statistical systems on random surfaces. 
The range of apphcation of these ideas include non-critical string theory as 
well as Quantum Chromodynamics (QCD) in the large limit. Progress in 
this area has been possible because the mathematical knowledge on random 
matrices has increased dramatically in the last fifteen years [0]. 

An important mathematical object that appears naturally in the discus- 
sion of random matrices is the integral over the unitary group 0. This 
integral has been applied to the solution of the Two matrix model 0, 
and, more recently, to the Migdal-Kazakov model of "induced QCD" 

On the other hand, there is considerable expectation that supersymmetry 
might play an important role in the physical world ||^. Since there is a natural 
extension of matrices to supermatrices |^ , it is important to understand the 
properties of random supermatrices and study their relevance for the theory 
of random surfaces. 

In this paper we will start the study of random supermatrices related to 
a supersymmetric extension of the Itzykson-Zuber integral. In particular we 
will calculate some non-vanishing correlations of superunitary matrices. The 
paper is organized as follows : section 2 contains a very brief review of some 
basic properties of Grassmaniann manifolds, which will be used in the sequel. 
Section 3 contains the main result of this paper, which is the calculation of 
the IZ integral for U{m\n). This result is subsequently used to find some 
particular correlations among the supermatrices elements. 
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2 Basic properties of a superspace 



Since the purpose of this paper is to extend the Itzykson-Zuber integral to the 
case of the supergroup f/(m|n), which can be described by matrices acting on 
a superspace (supermatrices), we briefly review some of the basic properties of 
the linear algebra together with the differential and integral calculus defined 
over a Grassmann algebra. This sets the stage for the next section and also 
fixes our notation. For a more detailed and complete reference on these 
matters the reader is referred to Ref. @]. 

Let us consider a superspace with coordinates = {q\9°'),i = 1, . . . m, 
a = 1, . . . ,n such that g* {6") are even (odd) elements of a Grassmann al- 
gebra. This means that z^z'^ = (^—lyi^^+'^C^) z^ , where e(-P) is the Grass- 
mann parity of the index P defined by e{i) = 0, mod{2); e{a) = 1, mod{2). 
Also we have that e = J2^{Pi)- The above multiplication rule 

implies in particular that any odd element of the Grassmann algebra has zero 
square, i.e. it is nilpotent. 

Supermatrices are arrays that act linearly on the supercoordinates leaving 
invariant the partition among even and odd coordinates. To be more specific, 
the supercoordinates can be thought as forming an (m + n) x 1 column vector 
with the first m entries (last n entries)being even (odd) elements of the 
Grassman algebra. In this way, an (m + n) x (m + n) supermatrix is an array 
written in the partitioned block form 



where the constituent matrices have the entries Aij, Bi^, CaijDap- Be- 
sides, Aij, Da/3 {Bia, Cai) are even (odd) elements of the Grasmann algebra 
in such a way that the parity array of the supercoordinate vector columnn is 




(1) 
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preserved. The parity of any supermatrix element is e(MpQ) = e(P) + e{Q) 
and defines the muhiphcation among supermatrix elements. The addition 
and multiplication of supermatrices according to the rules 

(Ml + M2)PQ = (Mi)pQ + (M2)pQ, {MiM2)pQ = {Mi)pr{M2)rQ, 

is such that it produces again a supermatrix. The inverse of a supermatrix 
can be constructed in block form, in complete analogy with the classical case 
and it is well defined provided and D^^ exist. The inverse of an even 
matrix is calculated in the standard way. 

The basic invariant of a supermatrix under similarity transformations is 
the sup ert race 

m+n 

Str{M) = Tr{A) - Tr{D) = ^ (-1)^(^)Mpp, 

p=i 

which is defined so that the cyclic property Str{MiM2) = Str{M2Mi) is 
fullfilled for arbitrary supermatrices Mi,M2. The above definition of the 
supertrace leads to the construction of the superdeterminant in the form 
Sdet{M) = exp[Str{lnM)], which is explicitly given in the following two 
equivalent forms 

Sdet(M) - ^^^(-4 - BD-'C) _ detjA) 
bdet{M)- ^^^^^^ - det{D-CA-^By 

The above expression is written only in terms of even matrices in such a way 
that the determinant has its usual meaning. The superdeterminat has the 
multiphcative property Sdet{MiM2) = Sdet{Mi)Sdet{M2). 

The definition of the adjoint supermatrix follows the usual steps by re- 
quirmg the identity {y^*MpQZ^y = z^* pqV^, for an arbitrary bilinear 
form in the complex supercoordinates y^, where * denotes complex conju- 
gation. Since the usual definition of complex conjugation in a Grassmann 
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algebra {y^y^Y = y^*y^* reverses the order of the factors without introduc- 
cing any sign factor , we have the result M^pq — Mqp* as in the standard 
case. 

A hermitian {m + n) x (m + n) supermatrix M is such that — M 
and it has (m + nY real independent components. The following proper- 
ties are also fuUfilled: (i) (Mt)t = M, (ii) (MiMs)^ = M2^mY and (iii) 
Sdet{M'^) = Sdet{MY- A unitary (m + n) x (m + n) supermatrix U is such 
that UW = WU = I { where / is the identity supermatrix) and also has 
(m + nY real independent components, which have the additional property 
that {S detU) {S detUY = 1- The set of all (m + n) x (m + n) unitary super- 
matrices form a group, called the supergroup U{m\n), under the operation 
of supermatrix multiplication. 

In the next section we will deal with derivation and integration with 
respect to the elements of a supermatrix, which will be considered as su- 
percoordinates there. To this end, we summarize here the basic properties 
involved, in terms of the supercoordinates z^. Differentiation and integration 
over the even supercoordinates follow the same rules and properties as the 
corresponding operations with complex numbers. On the other hand, the sit- 
uation with respect to the odd supercoordinates is greatly simplified because 
of the nilpotency property, which leads to the conclusion that an arbitrary 
superfunction F{q\ 9°') can be expanded into a finite set of products of the 
odd supercoordinates. In fact we have 

F{q\ 0^) = f{q') + + fa^a^Qy^'O'^- + ... + ■ ■ ■ (3) 

where all the functions fai...ak are completely antisymmetric in all the subindices. 
Also, a particular odd supercoordinate 9°''' can only appear linearly in the 
above expansion. Right and left derivatives with respect to the odd coordi- 
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nates are defined according to 

QRp QLp 

dF{q\e^) = ^de- = de-^, (4) 

where the operator d has zero Grassman parity and acts distributively upon a 
product of odd coordinates: c/(0"i0"2 ,,,0^k)= rfr^^^^ , , , Qa^^Qa^^QCi ^q^^^ 

QaiQa2 dO"'''. Both, right and left derivatives satisfy adequate Leibnitz rules 
which can be directly obtained from the definitions (^). Clearly, for even co- 
ordinates right and left derivatives coincide. 

Integration over odd Grassmann variables is defined, according to Berezin, 
by the basic properties ||^ 

J dO'^ = 0, j dO'^e'^ = s'^p, (5) 

which allows the calculation of integrals over arbitrary functions using the 
expansion @ together with the multiplication rules of supernumbers. The 
full integration measure of the superspace under consideration is given by 
[dz] = dq^ . . . dq^dO^ . . . dO"^, where a particular order in the odd differentials 
has been chosen. Under a change of supercoordinates z'^ = z'^{z^) the 
integration measure transforms as 



[dz'] = Sdet{-^)[dz]. (6) 

3 Integration over the unitary supergroup U(m\n 

In this section we calculate the extension to (m + n) x (m + n) supermatrices 
of the IZ integral. Our method, mutatis mutandis, follows closely to that of 
Itzykson and Zuber 0. Let us consider the integral 

/(Ml, M2; P)= f [rff/]e'^^*^'(^^^^*^^^'), (7) 



where Mi, M2 are hermitian supermatrices which can be diagonahzed and 
j3 is an even parameter. Up to a normahzation factor yU (to be fixed later), 
we define the integration measure over U{m\n) by 

[dU] = /i n dUpQdU*pQ 5{UU^ - I). (8) 

P,Q=1 

Here the 5-function really means the product of {m + unidimensional 5- 
functions corresponding to the independent constraints set by the condition 
UU'^ = I. In the case of an odd Grassmann variable 6, 

5(e -e) = {e-e) = j dne^'^^-^l 

where vr is another odd integration variable. It is important to observe that 
the measure in (0) possesses 2mn real independent odd differentials. 

The measure (j^) is invariant under independent right and left multipli- 
cation by arbitrary unitary supermatrices. In this way, / depends only upon 
the eigenvalues of Mi and M2, which are given in the corresponding diagonal 
supermatrices Ai and A2. Our notation is such that the first m eigenvalues 
of A are identified by Aj, while the remaining n eigenvalues are denoted by 
Aq. Such partition is characterized by the following parity assignment of the 
eigenvector components Vp,Vp : e(Vp) = e(P),e(Vp) = e(P) + 1. Before 
sketching the calculation of the basic integral given in Eq. (|^), we state our 
final result and make some general comments. We obtain 

/(Ai,A2;/3) = S(Ai,Ai)S(A2,A2) x x 

Here, A is the usual Vandermonde determinant 

A(A) = n(A.-A,), A(A) = n (A. - A^) (10) 

i>j a>f3 
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and the new function that appears is 



m n 



(11) 



i=l a=l 



We observe that the polynomial S(A, A) is completely symmetric under in- 
dependent permutations of A, A. 

The expression is completely determined up to a normalization factor 
related to that of the measure in equation (|^). This situation is analogous to 
the standard IZ case where the required factor can be fixed directly from the 
corresponding expression by taking the limit Ai, A2 — > in a convenient way 
and demanding J[dU] = 1, for example. This procedure leads to the correct 
factors in Eq.(3.4) of Ref. 0. In our case, a similar limiting procedure leads 
to the conclusion that J[dU] = 0, precisely due to the appearance of the 
S(A, A) functions in the numerator. This is not an unexpected result since 
we are dealing with odd Grassmann numbers. For this reason we have chosen 
the normalization factor in such a way that 



constructing the superunitary invariant Laplacian operator on hermitian su- 
permatrices, which is given by 



where the derivatives are both either left or right derivatives. This is a 
particular case of the superlaplacian constructed for curvilinear coordinates 
and corresponds to the metric ds'^ = Str{dM'^) which is invariant under 



/(Ai, A2; P) = S(Ai, Ai)S(A2, A2)/3'""/(Ai, A2; P)I{Xi, A2; 



where I{di, ^2; P) is the corresponding IZ integral. 

Now we give some details of the proof of our result (^). We begin by 




(12) 
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global unitary transformations dM UdMU'^ for hermitian supermatrices 
M. 

Next we consider the propagator 

f{Ml,M2;t)=<M^\e'i^'\M2>, (13) 
which, by construction, satisfies the differential equation 

[^t-lDi^fiMuM2;t)=0, (14) 
together with the boundary condition 

/(Mi,M2;0) =< M1IM2 >= 5(Mi - M2), (15) 
at t=0. It can be shown that the propagator / is given by 

/(Mi,M2;t) = T7-^e^^*''(^^-^^)', (16) 

in analogy with the standard case. Since the supertrace of a squared su- 
permatrix is not positive definite, the factor i of the exponential in (|13D is 
included to insure convergence in further manipulations. 

Let us now consider the time evolution of a wave function g{Mi; t) satisfy- 
ing ([T^) and such that at t = reduces to a known function g{Mi] 0) = g{Mi), 
which is invariant under U{m\n) transformations. As a consequence, ^ is a 
function of the eigenvalues A's of M only and it is symmetric under the sep- 
arate permutations of the sets {Aj} and {Aq,}. Such a wave function can be 
constructed using the propagator in Eq.(|l^), as 

~g{M,; t)= j [dM]f{M,, M; t)^(M; 0), (17) 

where the integration measure over hermitian supermatrices is given by 
[dM] = lip dMpplls ji^s dMpsdM^g and it is invariant under U{m\n) trans- 
formations. 
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Now, we make a change of integration variables, from the initial M's to 
radial (A) and angular (U) variables given by M = UAW. This change of 
variables leads to 

~g{M,-t) = J [dA][dU]Jf{M,,UAU^-trg{A;0), (18) 

where the jacobian J is to be determined. The jacobian factorizes in the 
form J = A^(A)Ji(t/). The piece Ji{U) will remain incorporated to [dU] 
while we will find explicitly. Let us observe that according to Ref.0 
there are some restrictions to diagonalize a hermitian supermatrix. We are 
assuming that such forbidden points constitute a set of zero measure in the 
configuration space defined by {Mpq}. In virtue of the invariance of / under 
U{m\n), with the particular choice arising from Mi = \^Ai\^''', together with 
the invariance of [dU] under right and left multiplications, we arrive at the 
result 

~g{A,-t) = I [dA][dU]A\A)fiA,,UAW;trgiA), (19) 

which shows that the wave function preserves its invariance under U{m\n) 
transformations. 

The next step is to construct the wave function 

e(A;t)=A(A)^(A;t) (20) 

and to realize that 

k{A,,A-t) = A(Ai)A(A) J [rff/]/(Ai,f/Af/t;t), (21) 

is the kernel that corresponds to the propagator of ^{A;t). 

Before calculating such kernel in an independent way, so that we can use 
Eq. (|2T| ) to find the value for the angular integration, we now proceed to the 
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calculation of the jacobian A^. This calculation is equivalent to finding the 
volume element in spherical coordinates. The general property to be used 
here originates in the differential geometry of supermanifolds and states 
that given a metric tensor such that the length invariant is ds^ = dz^ gpgdz'^ , 
the correct integration measure over the manifold is given by ^/g [dz] , where 
g = Sdet{gpQ). As we mentioned previously, the adequate length element in 
the case of hermitian supermatrices is ds^ = Str{dM^). After making the 
change of variables M = U AW we are left with 

ds^ = Str{dA^) + Str{[U+dU, A]^). (22) 

Since we are only interested in extracting the A dependent piece of ^ = 
yjQkQu we can further consider the representation U = e*^, in terms of 
another hermitian supermatrix H. Taking H to be infinitesimal in such a way 
that WdU = idH., we obtain 

m n 

ds^ = ^ 5ijdXidXj — ^ 5af3d\adXi3 + 

i,j=l Q!,/3=l 

m n 

+ ^{^a — hYdHabdHha + ^ (A^ — XisYdHaisdHpot — 

a,b=l a,(3=l 
n,m m,n 

~ X! (-^a ~ ^pfdHapdHpa — ^ (Xa " XbYdHabdHha, (23) 

a,/3=l a,b=l 

for the expression (|2^) . From here we conclude that 

. ^ ^ n^>J{K-XJ)Ua>f3{K-Xp) ^ A(A)A(A) 

^ ^ U^AK-K) S(A,A) • ^ ^ 

Having completely determined the wave function ^(A;t) of eq.(^), we 
now calculate its equation of motion. The starting point is eq.(|l4|) for the 
function ^(A; t) = ^/A. Also we will need the expression of the superlaplacian 
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in curvilinear coordinates, which is given by 



1 f)L aL^ 

71 ^J-'''^'^''^ 



(25) 



where g^^ is the inverse of the metric tensor gpQ. The expression (|25|) further 
simphfies in our case, since we are applying the superlaplacian to a function 
which does not depend on the angular variables. In this way we obtain 



daA;t) I 
dt 2A 



tidK^ dK^A^' hdXa 9A/a^' 



(26) 



for the time evolution of .^(A; t), where we have dropped the refference to left 
derivatives since all variables are even. The above equation reduces to 



dt 2 



d' 



f(A;t), 



(27) 



in virtue of the following property of the function A 



A = 0. 



(28) 



The above equation can be proved from the expression 

1 



7^1 d\ 



+ 2E 



1 



(Aj - Xa){K - Xj) (Aj - Xa)' 



+ 



i,a,l3^a (Ai - Aa)(Ai - Xp) 



(29) 



together with the analogous relation 

1 



0=1 ^-^a a,i 



(Aj — Aq,)^ 



+ E 



(Aj — AQ,)(Aj — Aq.) 



+ 



+2 — - — —] 

a,P^a,i (Ai - A„)(A„ - A/3) 



(30) 
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The calculation of the above equations (^) and (^) has already used the 
identity 

1 1 1 



+ 



0, 



(/il - yU2)(Ail - /is) (/i2 - /i3)(/il - At2) (/is - Ail)(/i3 - /i2) 

for the cases /i = A, A. The difference in the LHS of Eq. f^SD further reduces 
to 

/ m q2 n Q2 



XI Q\ 2 XI 



A = A 



with 



A. 



Aa + Afl — 2Aj 



o/3j 



5i 



(At — Aa)(Ai — A/3)(Aq — Xp) 



Aj + Aj — 2Aq, 



(31) 



(32) 



(33) 



(Aj — Aa)(Aj — Aj)(Aj — Aq.) 
The result given in Eq. (|28|) is finally obtained due to the antisymmetry prop- 



erties An 



The last step in the proof of Eq.(^) is the observation that the kernel 
K{Ai, A2]t) in eq.(|2T|) satisfies the following differential equation 

Wf(Ai.A,;t) ! a= V «M K(^ A ri41 



with the initial condition 



ir(Ai,A2;0) = (5(Ai-A2). 



(35) 



Moreover we require that K{Ai, A2;t) be separately antisymmetric under 
permutations of the {Aj} and the {Aq,}. Then, the solution of Eq. (|3^) is 

1 11 



i^(Ai,A2;t)) 



(27rit)(™+")/2 m! ra! 



(iet(e 



j(Ai,i-A2j)2/2t 



)det{e 



-i{Xi,c-X2,0r/'2t 



), 

(36) 



12 



which has the same symmetry properties of ^(A;t). From Eqs.(|T6|), (|2T| ) and 
we obtain 



I df/e-^^*^(^^^^^^') = factor x ("i 



m(m — l) n(n — l) 
mn / ^ \ ^-71 / ^ \ ^^"2 — - 



1\ 1\ 2 / 1 
[ul [~ it 



A(Ai)A(A2) ' ^ ^ 



where factor is independent of Ai and A2. Now, we use that e'^*'^^'"^'^^) 
Sdet{A) with A = e-^^'^^+^a) to get 



m(m — 1) , ^ nin — V) 



X , I ) 

A(Ai)A(A2) 

where (3 = j^. Factor is an arbitrary constant since it depends on (see 
eq.(|^)). We fix it in such a way that equation (9) holds. The imphcations of 
our formula (P) for the character expansion on U{m\n), together with a de- 
tailed calculation of the correlators of superunitary matrices will be discussed 
elsewhere. Here we content ourselves with the following result 

jm [Str(KmiU^)f = { ° 1; ^ : ■ • - 1 (39) 

The simplest identities that follow from the above result are 

J[dU]Ui,j,Ul,^...U,^,^Ul,^=0, forr = l,2...mn-l, (40) 

[[dU]UuUl...UmmUl^ = ^, forf/e[/(m|l), (41) 
J ml 

/ [dU]U22Ul2...U^+n 1+nUUn 1+n = for U G U {l\n) . (42) 
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